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ABSTRACT 

We use three-dimensional hydro-dynamic simulations to investigate the effects of a transient pho- 
toionizing ultraviolet (UV) flux on the collapse and cooling of pregalactic clouds. These clouds have 
masses in the range 10 5 - 10 7 Mq, form at high redshifts (z > 18), and are assumed to lie within 
the short-lived cosmological HII regions around the first generation of stars. In addition, we study 
the combined effects of this transient UV flux and a persistent Lyman- Werner (LW) background (at 
photon energies below 13.6eV) from distant sources. In the absence of a LW background, we find that 
a critical specific intensity of Juv ~0.1x 10~ 21 ergs s _1 cm~ 2 Hz -1 sr _1 demarcates a transition from 
net negative to positive feedback for the halo population. A weaker UV flux stimulates subsequent 
star formation inside the fossil HII regions, by enhancing the H2 molecule abundance. A stronger UV 
flux significantly delays star-formation by reducing the gas density, and increasing the cooling time, 
at the centers of collapsing halos. At a fixed Jtjv, the sign of the feedback also depends strongly 
on the density of the gas at the time of UV illumination. Regardless of the whether the feedback 
is positive or negative, we find that once the UV flux is turned off, its impact stars to diminish af- 
ter ~ 30% of the Hubble time. In the more realistic case when a LW background is present, with 
Jlw > 0.01 x 10 _21 ergs s _1 cm -2 Hz -1 sr _1 , strong suppression persists down to the lowest redshift 
(z = 18) in our simulations. Finally, we find evidence that heating and photoevaporation by the 
transient UV flux renders the ~ 10 6 M Q halos inside fossil HII regions more vulnerable to subsequent 
H2 photo-dissociation by a LW background. 
Subject headings: cosmology: theory - early Universe - galaxies: high-redshift - evolution 



1. INTRODUCTION 

Semi-analytical and numerical studies agree that the 
first generation of stars are likely to have formed at very 
high redshifts, z > 20, at locations corresponding to rare 
peaks of the fluctuating primordial density field. Numer- 
ical simulations suggest that the first collapsed objects 
host high-mass (~ 100M Q ), low-metalicity (so-called 
'Population III', hereaft er 'PopIII') stars (Abel et al. 
l2002tlBromm e t al. 2002J). In the absence of any feedback 
processes, these stars, and their accreting remnant black 
holes (BHs), could significantly reionize the intergalactic 
medium (IGM). Recent evidence from the z ~ 6 quasars 
discovered in the Sloan Digital Sky Survey (SDSS), 
whose s pectra exhibit dark and sharp G unn-Peterson 
troughs i|Fan et all2006tlMesinger fc Haimanl2004 - and 
from the cosmic microwave polarization anisotropies 
measured by the WMAP satellite, which imply a conser- 
vative electron scattering optical dept h r e = 0.09 ± 0.03 
IjPaee et al . 2006; Spergel et al.l2006q) . suggest that cos- 
mic reionization was delayed to a later stage of structure 
formation, z ~ 6 - 10. Such a late reionization would 
require a substantial suppressi on of very hi gh redshift 
ionizing sources (e.g., [Haiman fe Brvanl l^OOG 1 ). 

The impact that the first generation of stars 
would have on their surroundings plays a crucial 
role in how reionization progr esses a t high redshifts 
JHaiman fc Holderl l2003t ICenl 120031 IWvithe fc Loebl 
2003), but is poorly understood from first principles. 
Several feedback mechanisms are expected to be poten- 
tially important, including both direct chemical and en- 
ergetic input from the first supernovae, and radiative 
processes due to UV and X-ray radiation from the first 



stars themselves. The enhanced metallicity due to the 
supernovae will result in an increase in the cooling rate, 
leading to m ore efficient star formation and lower stellar 
masses l'e.g.. lOmukaill200rllBromm fc Loebl l2003). The 
signature of the metals that are produced may be seen in 
quasar absorption spectra, although semi-analytic work 
on the subject suggests that chemical feedback from pre- 
galactic objects did not play a large role in setting the 
observed intergalactic metallicity distribution at z < 6 
IjScannapieco et all 2002, 2003). In this paper we focus 
on radiative feedback, and postpone the study of PopIII 
metal pollution and feedback to future work. 

Radiative feedback can be either positive and nega- 
tive, in that it can enhance or suppress subsequent star- 
formation. Positive feedback can result when the en- 
hanced free-electron fraction from ionizing photons or 
from shocks catalyzes the formation of molecular hy- 
drogen (H2), which can provide the dominant cooling 
channel at high-densities and low temperatures. The 
catalyst electrons can be produced by X-rays emit- 
ted as a result of gas accretion onto early BHs (e.g., 
iHaiman et al.lll996j) . or by a previo us epoch of photoion- 
ization inside "fossil" HII regions llRicotti et alJ [2002b; 
IQh fc Haimanl2003HO'Shea et alJ2005lh or bv collisional 
ioniz ation in prot o galactic shocks (Shapiro & Kane 1987; 
iSusa et alJ 119981 lOh fc Haimar] |2002|) . Indeed, cos- 
mological simulations have noted net positive feedback 
close to the edge of H II regions (jRicotti et all 12002b: 
iKuhlen fc Madaul I2005D . Negative feedback can result 
from chemical or thermodynamical effects. UV pho- 
tons in the Lyman- Werner (LW) band of H2 can dis- 
sociate these molecules, there by reducing their effective- 
ness in cooling the gas (e.g., IHaiman et alJll997l 2000; 



iCiardi et al1l2000HMachacek et all 1200 11. Active radia- 
tive heating can photo-evaporate ga s in low — mass halos 
iEfstathioul 119921: JBarkana fe Loeblll999t iShapiro et all 
2004). Additionally, a past episode of photoioniza- 
tion heating in inactive "fossil" HII regions can leave 
the gas with tenacious excess entropy, reducing the gas 
densities, h indering H2 formation, cooling and collapse 
(IQh fc Haimanll20ml) . 

The above feedback effects, their relative importance, 
and their net outcome on star formation within the pop- 
ulation of early halos, is not well understood ab-initio, 
and is poorly constrained by observations at high red- 
shifts. Although invaluable in furthering our under- 
standing of t he main physical concepts, semi-analytic 
studies (e.g.. iHaiman et alJl!997t IQh fc Haimanl l2003t 
IMacIntvre et alJ I2005f) do not fully take into account 
the details of cosmological density structure and evolu- 
tion, which can be very important. On the other hand, 
numerical studies are limited in scope due to computa- 
tional restrictions associated with full radiative transfer 
on such large sca les and large dynamical ranges (see, e.g., 
Iliev et al. 2006 for a recent review). In particular, it is 
difficult to couple radiative transfer and hydrodynamics 
- most work to date has focused on e ither one or the 
other issue (with some exceptions, e.g.,|G ncdin fc Abell 
2001; Ricotti et al. 2002a: Shapiro et al. 2004). Tech- 
niques approximating full radiative transfer can provide 
a crucial speed-up of computing time, but such simula- 
tions still do not provide a large statistical sample for 
a detailed study of radiative feedback, and/or can be 
limited by a small dynamical range, thereby missing the 
smallest halos whic h would be most sus ceptible to nega- 
tive feedback (e.g.. lRicotti et al1l2002aj) . 

The purpose of the present paper is to statistically in- 
vestigate UV radiative feedback associated with the first 
generation of stars. Prior numerical studies have ei- 
ther focused on radiati on in a single differen t band, 
such as LW photons llMachacek et alJ 12001 ) or X- 



rays l|Machacek et al.ll2003t iKuhlen fc MadauN2005|) or 
lacked quantitative statistics and have not included 
photo-heat ing and photo -evapo ration of lo w-mass halos 
JRicotti et alJl2002alfi lO'Shea et al"ll2005|) . The recent 
work bv lAlvarez et alJ 1J2006J) has studied the impact of 
photoionizing radiation in detail within a single HII re- 
gion, but without self-consistently modeling the hydro- 
dynamics. In the present study, we quantify the com- 
bined effects of UV photo-ionization and LW radiation 
from nearby Pop III star formation. Rather than sim- 
ulating the radiative transfer within an individual HII 
region, we take a statistical approach. In particular, we 
examine how halos which are in the process of collapsing 
are affected by spatially constant but potentially short- 
lived radiation backgrounds at various intensities. This 
allows us to calibrate the sign and amplitude of the re- 
sulting feedback, in order to include these effects in future 
semi-analytic studies. To this end, we carry out simula- 
tions in which a large region is photo-ionized for a short 
period of time, and neglect radiative transfer effects (we 
discuss the impact of this approximation in more detail 
in §|3 below). 

The rest of this paper is organized as follows. In § [3 
we describe the simulations. In § |3 we present the re- 
sults of the simulations with photoionization heating, but 
without a LW background. In §0] we discuss simulation 



runs that also include a LW background. Finally, in §0 
we summarize our conclusions and discuss the implica- 
tions of this work. For completeness and for reference, 
in the Appendix, we present the dark matter halo mass 
functions found in our simulations. 

Throughout this paper, we adopt the background cos- 
mological parameters (51a, ^m, ^b, n, o~a, Hq) = (0.7, 
0.3, 0.047, 1, 0.92, 70 km s^ 1 Mpc -1 ), consistent with 
the measurements of the power spectrum of CMB tem- 
perature anisotrop ies by the first yea r of data from the 
WMAP satellite l|Spergel et all [20031. The three-year 
data from WMAP favors decrea sed small-scale powe r 
(i.e. lower values for as and n s ; iSpergel et al J 12006a). 
which would translate to a ~ 15% redshift delay, but 
would not change our conclusions. Unless stated other- 
wise, we quote all quantitities in comoving units. 

2. SIMULATIONS 

We use the Eulerian adaptive mesh refinement (AMR) 
code Enzo , which is described in greater detail elsewhere 
l)Brvanlll999t iNorman fc Brvanlll999j) . Our simulation 
volume is I (h^ 1 Mpc) 3 , initialized at z\ n \t = 99 with 
densit y perturbations drawn from the lEisenstein fc Hul 
(1999) power spectrum. We first run a low resolution 
(128 3 particles), dark matter (DM) only run down to 
z = 15, to find the highest density peak in the box. 
We then re-center the box around the spatial location 
of that peak, and rerun the simulations with the in- 
clusion of gas and at a higher resolution inside a 0.25 
h^ 1 Mpc cube, centered in the 1ft, -1 Mpc box. This re- 
fined central region has an average physical overdensity 
of <5(zi n it) = p/p — 1 = 0.1637, corresponding to a 2.4 a 
mass fluctuation of an equivalent spherical volume. We 
use such a biased region for our analysis since it hosts 
a large number of halos at high regshifts. This not only 
provides good number statistics, but also helps mimic 
a pristine, unpolluted region which is likely to host the 
first generation of stars. We postpone a lower redshift 
analysis to a future work. 

Our fiducial runs are shown in Tabled Our root grid 
is 128 3 . We have two additional static levels of refine- 
ment inside the central 0.25 h~ x Mpc cube. Further- 
more, grid cells inside the central region are allowed to 
dynamically refine so that the Jeans length is resolved 
by at least 4 grid zones and no grid cell contains more 
than 4 times the initial gas mass element. Each addi- 
tional grid level refines the mesh length of the parent 
grid cell by a factor of 2. We allow for a maximum of 
10 levels of refinement inside the refined central region, 
granting us a spatial resolution of 7.63 h~ x pc. This co- 
moving resolution translates to 0.36 h^ 1 proper pc at 
z = 20. We find that this resolution is sufficient to 
adequately resolve the gross physical processes in the 
few xlO 5 — I0 7 Mq halos of interest in this work by 
comparing with higher mass and spatial resolution runs 
(not shown here). The dark matter particle mass is 747 
Mq. We also include the non-equilibrium reaction net- 
work of nine chemical species (H, H + , He, He + , He ++ , 
e~, H 2, H^, H~) using the algorithm of lAnninos et alJ 
(1997), and initialized with post-recombination abun- 
dances from Anninos & Norman (1996). Our analysis 
below is based on the central refined region; the low 
resolution dark matter outside the refined region serves 
to provide the necessary tidal forces to our refined re- 



gion. Readers interested in further details concerning 
the simulation metho dolog y are encouraged to consult, 
e.g. JMachacek et all l)2001|K 

As shown in Table ^ we have performed five differ- 
ent runs without a LW background, distinguished by the 
duration or amplitude of the assumed UV background 
radiation (hereafter UVB), and six additional runs that 
include an additional constant LW background. For the 
UV radiation we assume an isotropic background flux 
with a T = 2 x 10 4 K blackbody spectral shape, normal- 
ized at the hydrogen ionization frequency, hvn = 13.6 
eV. The values of Juv are shown in Table H in units of 
10 _21 ergs s _1 cm -2 Hz" 1 sr _1 . The NoUVB run con- 
tains no UV radiation, and serves mainly as a reference 
run. The Heat0.08 and Heat0.8 runs include a UVB with 
Juv = 0.08, 0.8, respectively. The value of Juv = 0.08 
was chosen to correspond to the mean UV flux expected 
insid e a typical HII regio n surrounding a primordial star 
('e.g.. lAlvarez et aLll2006lh As we do not include dynami- 
cally expanding Hff regions in our code, the Heat0.8 and 
Flash runs can be viewed as extremes, corresponding to 
conditions close to the center and close to the edge of 
the HII region, respectively. More generally, studying 
a range of values of Juv is useful, since the UV flux 
of a massive, Pop-Ill star is uncertain. In the latter 
two runs, the UVB is turned on at zuvB.on = 25 and 
turned off at zuvB,off = 24.64. This redshift range corre- 
sponds to a typical theoretical stellar lif etime, ~ 3Myr , 
of a ~ 100M Q primordial (Pop-Ill) star l|Schaererll2002ft . 
The flash ionization run, Flash, instantaneously sets the 
gas temperature to T=15000 K and the hydrogen neu- 
tral fraction to shi = 10 -3 throughout the simulation 
volume, but involves no heating thereafter. This allows 
us to compare our results to those o flO'Shea et al. (2005), 
and to identify the importance of including the dynami- 
cal effects of the photo-heating. We also include an early 
UVB run, EarlyHeat0.8, with Juv = 0.8, zuvB,on = 33, 
and zuvB,off = 33.23, in order to study how our results 
vary with redshift (or equivalently, with the ionizing ef- 
ficiency of the first sources). Finally, the six runs in 
the bottom half of the Table repeat pairs of the NoUVB 
and Heat0.8 runs with three different constant LW back- 
grounds (Jlw = 0.001, 0.01, and 0.1, normalized at 
12.87eV in units of 10 _21 ergs s _1 cm~ 2 Hz -1 sr _1 , and 
assumed to be frequency-independent within the narrow 
LWband). 

We stress that we do not attempt to model reionization 
in this work. Rather, we focus on the statistical analysis 
of the feedback associated with the UV and LW back- 
grounds. In particular, we want to simulate the effect of 
a short-lived UV and a persistent LW background on a 
range of halo masses, in order to calibrate the net impact 
of fossil HII regions on subsequent star formation within 
these regions. In this we are aided by the large number 
of halos (few hundred) in our refined region, dozens of 
which manage to host cold, dense (CD) gas by the end 
of our simulation runs at z = 1 8. 

We use the HOP algorithm ifEiscnstci n & HutHl998l 
on the DM particles to identify DM halos. We then 
convert the resulting DM halo mass, Mom, to a total 
halo mass using the average conversion factor, Mhaio = 
Mum ^m/(^m — ^b)- We find that the halo masses 
defined in this manner agree to within a factor of two 
with masses obtained by integrating the densities over a 



TABLE 1 
Summary of Simulation Runs 



Run Name 


Juv 


ZUVB.on 


^UVB.off 


JvN 


Runs w 


ithout 


a LW background 




NoUVB 





NA 


NA 





Flash 


- 


25 


25 





Heat0.08 


0.08 


25 


24.62 





Heat0.8 


0.8 


25 


24.62 





EarlyHeat0.8 


0.8 


33 


32.23 





Runs with a 


LW background 




NoUVB 





NA 


NA 


0.001 


Heat0.8 


0.8 


25 


24.62 


0.001 


NoUVB 





NA 


NA 


0.01 


Heat0.8 


0.8 


25 


24.62 


0.01 


NoUVB 





NA 


NA 


0.1 


HeatO.8 


0.8 


25 


24.62 


0.1 



sphere whose radius is the halo's virial radius. 

In the analysis below, it will be useful to define the 
fraction of total gas within the virial radius which is cold 
and dense (CD), / c d. By cold, we mean gas whose tem- 
perature is < 0.5T v ; r , where r v ; r is the halo's virial tem- 
perature (for how virial tem peratures are associated with 
halos in the simulation, see Machacek et al. 2001). By 
dense, we mean gas whose density is > 10 19 M Q Mpc -3 
~ 330 cm -3 , roughly corresponding to the density at 
which the baryons become important to the gravitation 
potential at the core, taken to be an imme diate precursor 
to primordial star formation IjAbel et al, 2002). Hence- 
forth, we treat / c d as a proxy for the fraction of the halo's 
gas which is available for star formation. 

Additionally, we discount halos which have been sub- 
stantially contaminated by the large (low-resolution) DM 
particles outside of our refined region. Specifically, we re- 
move fromS our analysis halos with an average DM par- 
ticle mass greater than 115% of the refined region's DM 
mass resolution, 747M Q . Another possible source of con- 
tamination arises from closely separated halos. If some 
CD gas belonging to a halo is within another halo's virial 
radius (most likely in the process of merging) , the other 
halo could undeservedly be flagged as containing CD gas 
as well. To counteract this, we set / c d = for low-mass 
halos (< 2 x 10 5 M Q ) whose centers are less than ~ 5 h^ 1 
kpc away from the center of a halo containing CD gas. 

3. RESULTS WITHOUT A LW BACKGROUND 

In Figure ^ we show grey scale temperature projec- 
tions of a 20 ft. -1 kpc comoving region surrounding two 
halos nested in a filament. The rows correspond to the 
NoUVB, Flash, Heat0.08, and Heat0.8 runs (from top 
to bottom). The columns correspond to redshifts z — 
24.62, 24, 22, and 20 (from left to right). The scale is 
logarithmic, with black corresponding to T < 100K and 
white corresponding to T = 10 4 K. The halo in the lower 
(upper) part of each figure grows from M = 7.2 x I0 5 M© 
(M = 6.9 x 10 5 M Q ) at z = 24.62 to M = 1.4 x 10 6 M Q 
(M = 1.2 x 10 6 M Q ) at Z = 20, as measured in the 
NoUVB run. 

As one would expect, when the UVB is turned on, 
the gas is quickly ionized and heated. Gas which was 
previously at or close to hydrostatic equilibrium now has 
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Fig. 1. — Temperature projections of a 20 /i -1 kpc comoving region surrounding two halos nested in a filament. The rows, from top 
to bottom, correspond to the NoUVB, Flash, Heat0.08, and Heat0.8 runs. The columns, from left to right, correspond to redshifts of z = 
24.62, 24, 22, and 20. The scale is logarithmic, with black corresponding to T < 100K and white corresponding to T = 10 4 K. 



a greatly increased pressure gradient due to the increase 
in temperature. As a result, an outward-moving shock 
is formed, as clearly seen in Figure ^ f° r the last two 
rows (i.e. the runs which include a UVB with dynamical 
heating). Note that this shock is nearly absent in the 
Flash run. Subsequently, the gas in the dense filaments 
inside the shock is able to cool through Compton and H2 
cooling, and the shock stalls. The gas surrounding the 
halo starts infalling again. We explore these processes 
more quantitatively in § 13.11 

Note that the cores of the halos retain CD gas in all of 
the runs. In runs containing a UVB, the low density IGM 



outside of the filament still hasn't cooled below T ~ 10 3 K 
by z — 20; however, the filament itself shows evidence of 
positive feedback in the Flash and Heat0.08 runs, with 
lower temperatures at z < 24 than in the NoUVB run. 
Furthermore, it is evident that once the UVB is turned 
off, the dense filament is able to cool very rapidly, from 
T ~ 10 4 K to T ~ 10 3 K in Az < 0.6, due to the increased 
electron fraction, as we shall see below. 



3.1. Halo Profiles 



z = 24.62 



z = 23.00 




z=18.00 



10 10 2 

r (pc) 



10 10 2 

r (pc) 

= 24.62 



10 10 2 

r (pc) 



10 10 2 

r (pc) 

= 23.00 




10 10 2 

r (pc) 




10 10 2 

r (pc) 



10 10 2 

r (pc) 



Fig. 2. — Spherically averaged radial profiles of the same halo in the NoUVB (solid lines) and Heat0.8 (dashed lines) simulation runs. 
This halo was able to first form CD gas at z = 21 (z = 18) in the NoUVB (Hcat0.8) run, respectively. The left pair of figures show a 
snapshot at z = zrjVB.off = 24.62, the middle pair at z = 23, and the right pair at z = 18. The virial radius of the halo increases from 
-Rvir ~ 100 pc at z = 24.62, to i? v ir ~ 200 pc at z = 18. All quantities are shown in proper (not comoving) units. The Upper panels show 
the hydrogen density, mass- weighted gas temperature, gas cooling time, and radial velocity (clockwise from upper left). The Bottom panels 
show mass fractions of HI, HII, H2, and the number fraction of e~ (clockwise from upper left). 



To get a more quantitative idea of the feedback intro- 
duced by a UVB, in FigureElwe plot spherically averaged 
radial profiles for the same individual halo at redshifts 
z — 24.62, 23, and 18 (left to right), and in two different 
runs: NoUVB and Heat0.8 (solid and dashed lines, re- 
spectively). Figures in the top row show hydrogen den- 
sity, mass-weighted gas temperature, gas cooling time, 
and radial velocity (clockwise from upper left). Figures in 
the bottom row shows mass fractions of HI, HII, H2, and 
the number fraction of e~ (more precisely, f e is defined 
as the mass fraction of e _ , normalized such that each 
e~ is assumed to have the mass of hydrogen) (clockwise 
from upper left). The halo has a mass of M(z = 24.62) = 
3.42xl0 5 M Q andM(z = 18) = 2.37xlO 6 M [taken from 
the NoUVB run; note that the mass in the Heat0.8 run 
is somewhat smaller, e.g., M(z = 18) = 2.24 x 1O 6 M , 
due to photo-evaporation and a slight suppression of gas 
accretion as a result of the UVB] . 

From the profiles, one can see the impact of photo- 
evaporation in the Heat0.8 simulation run: the radially- 
outward moving shock mentioned above, as well as an 
accompanying decrease in density. As soon as the ioniz- 
ing radiation is turned off, the gas cools from a temper- 
ature of ~ 10 4 K to ~ 10 3 K quite rapidly, with the free 
electron number fraction (approximately corresponding 



to the bottom left panels of the bottom row of Fig. GJ , 
dropping two orders of magnitude by z = 23, Az = 1.62, 
since the UVB was turned off. Also, the shock starts 
to dissipate by z ~ 23, with most of the gas switching 
to the infall regime again. Despite the evident photo- 
evaporation, the presence of the UVB has caused over 
an order of magnitude increase in the H2 fraction, due to 
the increased (out-of-equilibrium) number of free elec- 
trons soon after zuvB,ori- 

We note that this halo managed to first form CD gas 
at z — 21 in the NoUVB case, but the formation of CD 
gas was delayed in the Heat0.8 case until z = 18. This 
delay can be crudely understood by looking at three fun- 
damental time-scales: the H2 cooling time, 
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Here, kn is the Boltzmann constant, T is the tempera- 
ture, Ah 2 is the H 2 cooling function, x e is the free elec- 
tron number fraction, n g , tihi, and uh 2 are the num- 
ber densities of all baryons and electrons, neutral hydro- 
gen, and H2, respectively, and 8 is the gas overdensity, 
5 = n g /n g - 1. 

The first column of Figure [21 shows a snapshot of our 
halos immediately prior to turning off the UVB. Note 
that the cooling time in the Heat0.8 run is several orders 
of magnitude lower than in the NoUVB run. 1 At large 
radii, this is due to Compton cooling, since the UVB 
dramatically increases x e and the Compton cooling time 
(cq. EJ) scales as x~ x . Additionally, the temperature 
increase to ~ 10 4 K allows for far more efficient line 
cooling from atomic hydrogen than in the NoUVB case. 
We also see that the structure of halos is important in 
accurately predicting feedback. Specifically, we note that 
the central region can behave quite differently than the 
outer regions of the halo. 2 

Initially, immediately after the radiation field is turned 
off, the gas cools due to atomic line cooling. However, 
this quickly becomes inefficient below temperature of 
about 6000 K. While the radiation is on, the H2 abun- 
dance is highly suppressed due to Lyman- Werner disso- 
ciation. However, as the temperature drops, the amount 
of H2 increases rapidly because of the high electron abun- 
dance. In fact the molecular hydrogen formation time is 
shorter than the recombination time and a large amount 
of molecular hydrogen is produced, xh 2 ~3x 10~ 3 , ir- 
respective ofjhe_densityand temperature (see the dis- 
cussion in Oh & Haiman 2002 for an explanation of this 
freeze-out value). 

Due to the relatively high densities throughout the halo 
(8 ~ 10 4 near the core), as well as the high initial x e , 
the majority of the hydrogen recombines shortly after 
zuvB,oft in the Heat0.8 case (i rC c(^uvB,off) < lMyr). Af- 
ter a few recombination times, n g ~ nni, and so eq. fT| 
can be simplified to 
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where xh 2 is the number fraction of H2. This highly 
efficient molecular cooling channel is largely responsible 
for quickly driving the temperature down to about 1000 
K (although a persistent LW background can counter this 
H2 -enhancement effect; see below). Hence, we note that 
the cooling times near the halo core are comparable for 
both runs at z = 23, with the cooling time in the Heat0.8 
case being larger by a factor of a few. 

This factor of only a few change in the cooling time is 
due to the remarkable cancellation of two strong effects, 
and can be understood by noting that the UVB induced 
photo-evaporation reduces n g by a factor of ~ 40 near 

1 The sharp drops in the cooling time in the NoUVB run cor- 
respond to annuli that include cold, low— density gas below the 
CMB temperature (~10 K), which is heated, rather than cooled, 
by Compton scattering. 

2 We do not include radiative transfer in our analysis, and the 
high-density regions (njj > 1cm -3 ), such as the cores of halos, 
migh t be able to self-shield against UV radiation lAlvarez et_al_J 
2006), decreasing feedback effects somewhat; see discussion in f 




Fig. 3. — Mass fractions of H2 at 2 = 18. Results are shown from 
the NoUVB {black crosses), Flash (blue dashes), Heat0.08 (green 
triangles), and Heat0.8 (red squares) simulation runs. 



the core at z — 23. Meanwhile, the UVB boosts xn 2 near 
the core by a factor of ~ 10. Since Compton cooling is 
ineffective at this stage, the cooling time is dominated 
by the H2 cooling channel, whose cooling time roughly 
scales as tu 2 oc (n g Xn 2 ) , given that the temperature 
is almost identical near the core at z ~ 23, and that the 
cooling function is very wea kly dependent on njj 2 in this 
regime (IQalli fc Pallalll99Sl) . 

From this crude estimate, one obtains an effective "de- 
lay" in the formation of CD gas in the Heat0.8 run with 
respect to the NoUVB run: 



/delay 
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40 
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This is in excellent agreement with the delay observed in 
the pair of simulation runs, where the halo obtains CD 
gas ~ 20 (60) Myr after our data point at z = 23 in the 
NoUVB (Heat0.8) case, yielding a delay of a factor of 

/delay ~ <J' 

3.2. Hi 'production 

As mentioned above, molecular hydrogen can provide a 
dominant cooling mechanism, especially in high density 
regions. To study the impact of the UVB on the forma- 
tion of H2 , we plot the H2 mass fractions as functions of 
halo mass in Figure |3| at the lowest redshift of our sim- 
ulations, z = 18. Results are shown from the NoUVB 
(black crosses), Flash (blue dashes), Heat0.08 (green tri- 
angles), and Heat0.8 (red squares) simulation runs. 

We note that in all of our runs which include a UVB, 
the molecular hydrogen fraction converges to a value 
which is nearly independent of the strength and dura- 
tion of the UVB. As seen in Figure |3J by the end of 
our simulations, most halos which have been exposed to 
a UVB in their past have settled on a value of a few 
x 10~ 3 for the mass fraction of H2. This freeze-out value 
is due to the fact that the number density of H2 follows 
its equilibrium value until about 3000 K, below which 
recombination proceeds more quickly than H2 formation 
and t he fraction freezes out at this point ( Oh & Haiman 
2002). This number is fairly independent of mass, though 
there is a weak trend towards higher mass fractions for 
higher mass halos, up to mass fractions of ~ 10~ 2 for 



M ~ 4 x 1O 6 M0. Note also that there is some evi- 
dence for a non-monotonic evolution of the H2 fraction, 
with mass fractions falling back down to ~ 4 x 10~ 3 by 
M ~ 3 x 10 7 M Q . 

In contrast, molecular hydrogen in halos which have 
not been exposed to a UVB (black crosses in Figure EU 
is distinctly sparser (over two orders of magnitude for 
M < 10 6 Af©) than in our other runs. Also, there is a 
stronger evolution with respect to mass, as well as more 
scatter (which makes sense, since the H2 abundance in 
this case is not a result of a freeze-out process, and de- 
pends strongly on local density and temperature). 

An interesting conclusion can be drawn from the sim- 
ilarity of H2 fractions in our runs which include a UVB. 
Namely, if our analysis of the dominant cooling processes 
in S 13. II is accurate, the differences between the CD gas 
fractions among our UVB runs is predominately due to 
disparate effectiveness of photo-evaporation. In other 
words, as the positive feedback (i.e. the increase of the 
xn 2 term in eq. 0J) is nearly independent of the strength 
and duration of the UVB in our runs, only the amount 
of negative feedback (decrease in n g ) causes variations in 
the delay in the formation of CD gas (see a more detailed 
discussion in §Q]below). We have also verified for several 
halos that this similarity in the total H2 fraction extends 
to the radial profiles of H2 . 



3.3. Ensemble Evolution of the Cold, Dense Gas 
Fractions 

As mentioned above, a quantity of particular impor- 
tance in studying the capacity of a halo at hosting stars 
is its cold, dense (CD) gas fraction, f c< \, defined above. 
Here we show the general trends for the evolution of this 
quantity for the population of halos in our simulations, 
as we vary the UVB. 

In Figure 01 we plot the total gas fractions (upper pan- 
els) and CD fractions (lower panels) as a function of total 
halo mass at redshift z = 18, the lowest redshift output of 
our simulations. The figures correspond to the NoUVB 
(top left). Flash (top right), Heat0.08 (bottom left), and 
Heat0.8 (bottom right) simulation runs. 

Note that while the total gas fractions of small halos 
(M < few xl0 5 M Q ) that have been exposed to a UVB 
are suppressed with respect to the NoUVB case, there 
is little immediate visual evidence of either negative or 
positive feedback for halos large enough to host CD gas 
(M > few xl0 5 M Q ). The Flash CD gas fractions show 
evidence of positive feedback in the mass range ~ 2 x 10 5 
- 10 6 Mq, while the Heat0.8 run shows evidence of strong 
negative feedback in the same mass range, with no ha- 
los hosting CD gas at M < 1O 6 M . These small halos 
are the ones which would be most affected by photo- 
evaporation. This lends further credibility to our asser- 
tion above that positive feedback in runs which include 
a UVB is fairly independent of the UVB strength, and 
hence the total feedback is set by photo-evaporation ef- 
fects (i.e. negative feedback). The Heat0.08 run has a 
near zero balance of positive and negative feedback. 

In order to better quantify the amount of suppression 
of CD gas in our models incorporating a UVB, as well 
as the evolution of such a suppression with redshift, we 
define the cumulative, fractional suppression of the halo 
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where V c N d oUVB (z) and N^ ni (z) are the total number 



f NoUVB / 

of halos with CD gas at redshift z in the NoUVB 
run and some given run i, respectively This expres- 
sion is well-defined for JV c N d oUVB (z) > iVf d oUVB (zuvB on); 
for V c N d ° UVB (z) = V c N d oUVB (z UV B,on), we set S N , cd (z) = 
0. Note that by definition, iV c N d oUVB (zuvB ; on) = 

-^cd m ( 2 UVB,on)- 

Similarly, we define the cumulative, fractional suppres- 
sion of the CD gas mass as 
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where M^ UVB (z) and M™ m (z) are the total mass of 
CD gas at redshift z in the NoUVB run and some given 
run i, respectively. The total CD gas mass is obtained by 
merely summing the CD gas masses for all of the halos 
in the simulation. As for equation JBJ, this expression is 
well-defined for M™ d oUYB (z) > M c N d oUVB (z UV B„„), and 
for M c N d ° UVB (z) = Af c N d ° UVB (z UVB ,on), we set S MtCd (z) = 
0. 

Equations JJJJ and (0 provide an estimate of how the 
CD gas has been affected by the presence of a UVB, 
following the turn-on redshift of the UVB, zuvB,on (the 
values at zuvb.oh are subtracted in order to provide a 
more sensitive measure of relative changes of CD gas). 
As defined above, <5jv,cd(z) = and 5m,c<±(z) — if the 
UVB has no effect. If the effect of a UVB is positive, re- 
sulting in positive feedback, 8n,ca(z) and 5m,ca(z) would 
be positive. If the effect of the UVB is negative, 5n,c&{z) 
and o~M,cd(z) would be negative. 

In FigureEl we plot the values of M^ m (z) (top panel) 
and N^ m (z) (bottom panel) in our four main simulation 
runs: rum = NoUVB (crosses), Flash (dashes), Heat0.08 
(triangles), and Heat0.8 (squares). The corresponding 
values of SM.cd(z) and <5jv, c d(z) are plotted in Figure H3 
in the top and bottom panels, respectively. The results 
are displayed for the Flash (dashes), Heat0.08 (triangles), 
and Heat0.8 (squares) simulation runs. Although some 
of the notable fractional changes shown in FigureH3might 
appear statistically insignificant due to the small number 
statistics inferred from Figure [3J it should be noted that 
these runs are not uncorrelated experiments. In other 
words, each of our runs in Tablets seeded with the same 
initial conditions, and so small relative changes compared 
to the NoUVB run are significant (i.e. the errors are not 
Poisson) . 

One can infer from Figures 03 and that the Heat0.8 
run shows evidence of strong negative feedback down to 
z ~ 20, with values approaching the NoUVB run by the 
end of our simulation (z = 18). Conversely, the Flash 
run exhibits strong positive feedback down to z ~ 20, 
and again approaches the NoUVB run by the end of our 
simulation. In the middle is the Heat0.08 run, which 
shows very little difference compared to the NoUVB run 
(initially there is some evidence of mild positive feedback 
down to a redshift of z — 21, but at redshifts below that, 
little evidence remains of a UVB ever being present). 

It is also interesting to note that while the halo num- 
ber in the Heat0.8 run shows negative feedback down 
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Fig. 4. — Total gas fractions (top panels) and cold, dense gas fractions (bottom panels) as a function of total halo mass at redshift z = 18, 
the lowest redshift output of our simulations. The four panels correspond to the NoUVB (top left), Flash (top right), Heat0.08 (bottom 
left), and HeatO.8 (bottom right) simulation runs. 



to z = 18 (bottom panels of Figures [3| and 0, the total 
mass of CD gas (top panels of Figures^ and\^ shows 
no such feedback at z < 19. The explanation for this 
apparent contradiction is that the total mass of CD gas 
is dominated by the largest halos (both because these 
halos are more massive and because the fraction of CD 
gas increases with mass), and as Figure 0] shows, these 
large halos are largely unaffected by the ionizing radia- 
tion. Conversely, the elimination of the CD gas from the 
lowest-mass halos even at z = 18 is a genuine effect (as 
is clearly visible in the lower right panel in FigQJ, but 
these halos do not contribute significantly to the total 
CD mass summed over all halos. 

Figure agrees well with the qualitative inferences 
drawn above. Furthermore, it explicitly shows that the 
critical UVB flux cutoff in our simulation between induc- 
ing a net negative and net positive feedback is Jtjv ~ 0.1. 
Halos which have been exposed to a fainter UVB ex- 
hibit positive feedback, whereas halos which have been 



exposed to a brighter UVB exhibit negative feedback. 
However, it is also important to note that any such feed- 
back is temporary, as all of our runs begin to converge 
by the end of our simulations at z = 18. The exception 
is that the HeatO.8 run shows persistent suppression of 
the smallest halos (with M < 10 6 M Q ) all the way down 
to z = 18. 

3.4. Relating initial densities at zuVB.on to subsequent 
suppression of cold, dense gas 

Here we attempt to generalize and physically motivate 
some of the results from the previous section. In par- 
ticular, we have already seen that feedback depends on 
Jtjv and A/haio- Here we examine whether a halo's ca- 
pacity for forming CD gas depends strongly on the prop- 
erties of its progenitor region at the time of the UV- 
illumination (ztjvb.oii)- Specifically, we expect those pro- 
genitor regions which are less dense at ztjvb.oii, and hence 
at an earlier evolutionary stage, to be more susceptible 
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Fig. 5.— Values of M™ ni (z) (top panel) and N™ l (z) (bottom 
panel) as defined in equations J7J and JHJ . The results are displayed 
for the NoUVB (crosses), Flash (dashes), Heat0.08 (triangles), and 
Heat0.8 (squares) simulation runs. 
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panel) as defined in equations J7J and JBJ. The results are de- 
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(triangles), and HeatO.8 (squares) simulation runs. 



to negative photo-heating and photo-evaporation feed- 
back than more dense regions. This is because the H 2 
photo-dissociation rate scales with the density, whereas 
H2-forming reaction rates scale with the square of the 
density; as a result, photo-dissociation becomes compar- 



atively more important at low densities l)Oh fc Haimanl 
2003). Below we focus on the HeatO.8 run, as it exhibits 
the strongest negative feedback. 

We divide the set of halos with CD gas at redshift z 
in our NoUVB run into two groups: those that also have 
CD gas in the HeatO.8 run (group 1), and those that 
do not also have CD gas in the HeatO.8 run (group 2). 
From Figure 0J one can note that at z = 18 it is possi- 
ble to define a rough mass scale that separates these two 
groups; namely halos with masses > 1O 6 M0 do not have 
their CD gas suppressed (group 1), and halos with masses 
< 10 6 M Q do have their CD gas suppressed (group 2). As 
stated above, we hypothesize that the physical distinc- 
tion between the two sets occurs due to their differences 
at the redshift they were exposed to the UVB, zuvb.oh- 
Specifically, we compare the mass- weighted, average den- 
sities of progenitor regions at zuvB.on = 25, which are to 
become our halos from groups 1 or 2 at some later z. We 
do this by tracing back all dark matter particles compris- 
ing each halo at z to their positions at z = zuvb.oii, and 
then obtaining the average gas density at that position. 

As there are too few halos to accurately construct the 
group 1 and group 2 mean density distribution functions 
(see bottom panel of Fig. |5J), we present their properties 
via a density cutoff. We adopt a simple criterion to de- 
fine a density cutoff, p C utoS, between the two groups. We 
chose Pcutoff so that the sum of the fraction of group (1) 
points below /5 C utoff and the fraction of group (2) points 
above p cu toS is minimized. Specifically, this fractional 
sum used as a proxy for the disjointness of the two dis- 
tributions is defined as 

/ fit = MIN [A(< p cutoff ) + / a (> Pcutoff)] (8) 

where /i(< p C utoff) is the fraction of halos in group 1 
which have mean densities less than the cutoff density, 
and fe(> /Ccutoff) is the fraction of halos in group 1 which 
have mean densities greater than the cutoff density. In 
essence, each term is the fraction of "misclassified" ha- 
los, and so we want to select p C utoff such that /fit is min- 
imized. The sum as defined in equation (JHJ) ranges from 

to 2, while our figure of merit ranges from /g t = for 
two completely disjoint distributions to /g t = 1 for the 
case where group 1 and group 2 are drawn from the same 
underlying distribution (not taking into account Poisson 
errors) . 

We plot values for the density cutoff (in units of the 
average comoving density, p) as a function of redshift 
in the top panel of Figure {7\ Our disjointness figure of 
merit is plotted in the bottom panel of Figure Note 
that for most redshifts, /fit <C 1, meaning that the group 

1 and group 2 density distributions are quite disjoint, 
and that the density cutoff, p C utoff, has a well-defined 
value. One can get a sense of the Poissonian errors 
associated with p C utoff by looking at the bottom panel 
of Figure since N group i ~ jv™^ - 8 - 10-20 and 
N grou p2 ~ AT c N°uvB_ iV Heato.8 _ 10 Qne should also note 
that /fit increases with decreasing p cu toff, which is to be 
expected as the density distributions can not have neg- 
ative values, so both distributions start being "packed" 
together as they approach zero. In other words, there 
is an intrinsic "noise" consisting of small environmental 
fluctuations (halo location, peculiar velocity, etc.), and 
this noise becomes more noticeable as p C utoff — > 0. In 
practice, it is difficult to disentangle this effect from an 
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Fig. 7. — The critical density Pcutoff, upper panel, for the pro- 
genitor gas at ZTjVB,on = 25 of halos that collapse at some later 
redshift z, roughly divides halos experiencing negative vs. positive 
feedback at redshift z (i.e. halos that were more/less dense than 
Pcutoff a t the time of illumination will experience positive/negative 
feedback) . p cu t ff is a mass- weighted mean density of the progen- 
itors pieces of the halo, shown in units of the average comoving 
density, p in our Heat0.8 run. The values of /g t , shown in the 
Lower panel, show a rough figure of merit for how well the fixed 
density p cu t ff separates halos into two disjoint categories (/fl t <K 1 
indicates a clear separation). See text and equation JSJ for defini- 
tions and discussion. 



actual merging of the two distributions. 

As expected, halos with less dense progenitor regions at 
zuvB,on are more susceptible to negative feedback. It is 
quite interesting to note that our density cutoff decreases 
exponentially with redshift, implying that an increasing 
fraction of the photo-heated mass will fall in the "bor- 
derline" region between negative and positive feedback. 
This provides further support for our earlier claim that 
the fossil HII region "forgets" the UVB as time passes. 
In other words, a strong UVB serves to merely delay 
the gas from cooling and collapsing; the gas eventually 
manages to cool, aided by an enhanced H2 fraction and 
enhanced infall (see Figure[3 and associated discussion). 
The length of this delay is a strong function of the den- 
sity of halo progenitor region s at zuvB.on, as one would 
expect from our analysis in § 13.11 

It is numerically impractical to run our simulations to 
rcdshifts much lower than z < 18, due to the rapidly in- 
creasing collapsed fraction in our refined region (see Fig. 
lAlfl . On the other hand, it would be interesting to know 
what eventually happens to most of the mass of our re- 
fined region. A step towards answering this question is to 
find out whether the majority of the mass of our refined 
region at z = ZuvB.on is located in overdensities below 
or above the lowest redshift p C utoff value shown in Figure 
With this motivation, we obtained a mass-weighted 
density distribution over randomly generated positions 
inside our refined region. We select a radius surrounding 
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Fig. 8. — Mass-weighted, cumulative density distributions for 
regions of M ~ 8.9 X 10 5 Af Q (solid curves) and M ~ 8.9 X 10 6 Af Q 
(dashed curves). Two redshift values are presented: ztjvb on = 33 
and 25, from left to right. Note that p is the comoving gas density. 



each position, such that a given number of DM particles 
lie within that radius (the number of particles is chosen 
to correspond to the halo masses capable of hosting CD 
gas). We then obtain a mass- weighted average density 
by averaging over the gas densities at the location of each 
DM particle inside our chosen radius. 

In Figure we plot the cumulative density distribu- 
tions (fraction of regions with mass- weighted density less 
than p) thus generated at z = 33 and 25, from left to 
right, and for regions of mass scale M ~ 8.9 x lO 5 Af0 
(~ 10 3 DM particles) and M ~ 8.9 x 10 6 M Q (~ 10 4 DM 
particles) with solid and dashed lines, respectively. Un- 
derstandably, the larger mass scales shift the mean den- 
sity towards larger values, due to the increased likelihood 
of averaging over dense patches. Also, we see that for 
regions of equal mass scales, the higher redshift counter- 
parts have a lower mean density, partly due to a smaller 
clumping factor and partly due to to the fact that we 
plot comoving density, which increases with decreasing 
redshift. 

Figure |H1 shows that the majority of the mass (> 90%) 
of our refined region is located in regions with mean den- 
sities lower than ~ 10/5, the lowest cutoff obtained by 
our analysis (see Figured). Hence, we cannot rule out 
the possibility of significant negative feedback at lower 
rcdshifts, not probed by our simulation. Nevertheless, 
we regard this as unlikely, for two reasons. First, ha- 
los will be centered around overdensities, not random 
points, and subsequent growth of the halo's mass need 
not be spherically symmetric; these effects will bias the 
relevant density distribution to higher values than shown 
in Figure HO Second, it is likely that most halos massive 
enough to host CD gas, which form in our biased re- 
gion at z < 18, already had a dense progenitor core at 
z = ztjvb. on- Indeed, we find that all halos hosting CD 
gas at z = 18 in the NoUVB run had some dense progen- 
itor gas (p > lOOp) at ztjvb, on = 25. Subsequent growth 
could be dominated by gas accretion onto these dense re- 
gions, and merging with other halos (rather than forming 
fresh halos entirely from low-density gas). Nevertheless, 
we emphasize that we chose to simulate a biased volume 
(favoring early collapse), and the above arguments will 
be less valid in a more typical patch of the IGM (with 



11 



10 7 - 



S io 6 



10 5 E" 



30 



20 



10 



: , , | , , , 










i 1 i i i = 

x NoUVB 


* X 










. EarlyHeat0.8" 


— ■ M 

E i x 








E 


2 X 








- 










.. 












■ X 








— 


: 


u 


x 




: 


_ 




. 


X 


_ 


: 








x : 






















_ 


X 










* NoUVB 












■ EarlyHeat0.8_ 




_ 












X X 








_ 


" 










XX 








- 










_ 












■ 










- 




v 




- 


" 








~ 


" 






X 


X " 






■ 






1 1 1 1 1 1 


1 




f 


— i — i — t — i — i — k — 



1 




_ I I 1 I I I 












- 


■ 












" 


■ 












- 


-1 






■ 




■ 


■ 




1 






























" 


■ ■ ■ ■ ■ ■ 












: 




- 




■ 


■ 






_ 


-1 


- 1 1 1 1 1 1 


1 




1 


■ 


■ 

1 


- 



20 



25 



30 



20 



25 



30 



Fig. 9.— Values ofM™ M (z) (top panel) and AT™ ni ( 2 ) (bottom 
panel) as defined in equations J7J and JBJ. The results are displayed 
for the NoUVB (crosses) and EarlyHeatO.8 (squares) simulation 



Fig. 10. — Values of Sm cd( z ) (top panel) and Sn c ^(z) (bot- 
tom panel) as defined in equations J3 and J^J, shown here in the 
EarlyHeat0.8 run. 



lower densities). 



3.5. Early UVB Run 



Ideally, one would want to explore all of parameter 
space by varying Jtjv, zuvb.oii and ZuvB.off, with differ- 
ent realizations of the density field. Unfortunately, given 
computational limitations, this is impossible. However, 
in order to confirm the trends we present above, we run 
another simulation, EarlyHeat0.8, in which we turn on a 
UVB, with an amplitude of Jtjv = 0.8, at zuvB.on = 33, 
and turn it off at ztjvb.oS = 32.23. We then repeat the 
analysis in § 13.41 The corresponding figures, Figures El 
and 1101 are presented below. 

Figure^]|shows that we once again find strong negative 
feedback down to z ~ 23. For z < 23, we see virtually 
no evidence of any feedback, lending further credibility 
to our interpretation above, that our other runs "forget" 
the episode of UV heating, and start converging to the 
NoUVB run by the end of our simulations (z = 18). 

I n Fig ure II 1 | we plot the density cutoff, ^cutoff , defined 
in § 13.41 for both the Heat0.8 (crosses) and EarlyHeat0.8 
(triangles) runs. For the sake of a direct comparison, 
this time we use physical units both for p C utoff, (proper 
cm -3 ), and for the time elapsed since the UVB turn-off 
(Myr) . Unfortunately, the drawback to having a simula- 
tion run with such an early heating episode is that there 
are fewer halos to analyze at earlier epochs. Specifically, 
in the epoch with evident negative feedback (z > 23), 
as seen below, there are only three redshift outputs con- 
taining both halos exhibiting suppression and halos not 
exhibiting suppression of C D gas (groups 2 and 1, re- 
spectively, defined in S I3.4JI . While it is difficult to draw 
strong conclusions from Figure ITT1 the density cutoff val- 
ues do appear similar in the two runs. 

We examined the radial profiles of the same halo pic- 



tured in Figure |2 to verify that we can apply the same 
cooling arguments as discussed in § 13.11 We compare 
the NoUVB and EarlyHeat0.8 runs at z = 30, shortly 
after zijVB.off = 33. As in Figure this redshift corre- 
sponds roughly to the regime where the induced shock 
begins to dissipate, and the gas starts falling back into 
the core. In all of the runs, the temperature drops to 
T ~ T v j r ~ 10 3 K near the core very soon after zuvB.off- 
As in S I3.1I we characterize the delay in formation of CD 
gas with (c.f. eq. |J5J) 
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From our simple cooling argument, we predict a nearly 
negligible delay for this halo. Indeed, the halo ends up 
forming CD gas at z = 20.5 in the EarlyHeat0.8 run, and 
at z = 21 in the NoUVBrun, showing a very small delay. 
Despite the halo's exposure to the UVB earlier in its 
evolution and subsequent lower gas density, the total 
negative feedback is reduced compared to the Heat0.8 
run. Compared to the Heat0.8 run, the negative feed- 
back, when expressed as the photo-evaporation term 
in the above equation ( n NoUVB / r jHcato.8-j j s sma ri er by 

a factor of ~ 2, and the positive feedback, when ex- 
presses as the H 2 fraction term in the above equation 
^Hoato.s^NoUVB) is larger by & factor of _ 2 . These 

changes are explained by more efficient Compton cool- 
ing (which more effectively eliminates the impact of the 
photo heating), and the lower gas density (which leads 
to a lower value for the H2 fraction in the NoUVB run, 
and hence a larger relative enhancemen t in the Heat0.8 
run), respectively IjOh fc Haimanl l2003) . 

3.6. The impact of not including radiative transfer 
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Fig. 11. — Values of p cu t ff (in proper cm 3 ) as a function of 
time elapsed since ztjvb off- Crosses correspond to our Heat0.8 run 
(i.e. ztjvb on = 25); triangles correspond to our EarlyHeat0.8 run 
(i.e. ztjvb, on = 33). 



Our simulations treat photo-ionization in the optically 
thin limit and so do not include radiative transfer ef- 
fects. This results in two differences compared to a self- 
consistent treatment. 

The most obvious effect is that all of our halos are 
ionized simultaneously, while in reality halos are ionized 
by very nearby stars with distances less than the few 
kpc radius of typical HII regions IjWhalen et al, 2004; 
iKitavama et al]l2004D . Nevertheless, we argue that the 
primary effect of this is to vary the flux felt by the halo 
and we explore a range of reasonable fluxes in our simu- 
lations. The exception is if the halo is so close that it is 
enveloped within the shock generated by the gas expelled 
from the halo hosting the star that produces the ionizing 
radiation. However, typically this shocked region occu- 
pies a volume of les s than 1% of the ionized region (e.g., 
Whalen ~alJl2004j) . 

The second, and more important, effect of radiative 
transfer will be to shield the high density cores of our 
minihalos. If the cores are not ionized, then both the pos- 
itive and negative feedback effect will clearly not occur 
in the neutral gas. lAlvarez et a l. (2006) estimate that 
self-shielding will set in at densities around a few par- 
ticles per cm -3 (depending on the strength of the flux 
and the size of the halo). This value is approximately 
the density we find in the cores of our simulated halos 
(e.g., Figure 2), and so we conclude that radiative trans- 
fer effects may play an important role in the cores of 
our halos. We note that at these densities, we typically 
find very little negative feedback anyway because of the 
short cooling times in the ionized gas. Most of the nega- 
tive feedback we observe arises due to the photo-heating 
of low-density gas which is then later accreted onto ha- 
los. This means that we do not expect our results to 
be strongly affected by the missing radiative transfer ef- 
fects. The effect, where important, will be to decrease the 
amount of feedback, making our statements about feed- 
back upper limits on the amplitude of the expected feed- 
back. Finally we note that, given time, the halos will be 
evaporated and eventually ionized despite the high den- 
sities in the core; however, this photo-evaporatio n time 
will be typically longer than 3 Myr (Haiman et al. 2001; 
IShaniro et aUEJOJlTliev et aTlEoO^ . 
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Fig. 12. — The figure illustrates the impact of the UV heating 
in the four different runs (Flash, Heat0.08, Heat0.8, EarlyHeatO. 8, 
as labeled in each panel). For each halo, we show, at z = 23 (or 
z = 30, in the early heating run), the ratio of the H2 fractions (red 
filled squares), the temperatures (blue empty triangle) and of the 
average gas density within the central 15pc (red filled circles) in 
the runs with UV heating, compared to the runs with no UVB. 
Note that the H2 fraction tends to increase by the same factor, 
regardless of the nature of the heating. As a result, the sign of the 
overall feedback (negative or positive) is determined primarily by 
the changes in the gas density, which does depend strongly on the 
type and amount of heating. 



4. THE IMPACT OF A LYMAN WERNER BACKGROUND 

While up to now we have ignored a possible Lyman- 
Werner background, such a background is likely to be 
present early on, and could have a strong impact on the 
H2 chemistry and gas cooling. In particular, the IGM 
is nearly optically thin, or quickly becomes so, at fre- 
quencies below 13.6eV (Haiman et al. 2000). For ref- 
erence, we note that one photon per hydrogen atom 
(the minimum UV background required for reionizing 
the IGM) would translate to a background intensity of 
J LW ~ 20[(1 + z)/21] 3 . Background levels 2-4 orders of 
magnitude below this value will be established well before 
reionization, and have the potential to alr eady photodis- 
sociate H2 molecules at thes e early epochs (Haima n et alJ 
H997tlMachacek et aill200H) . Furt hermore, and of more 
direct interest in the present study. lOh &: Haiman! (2003) 
have argued that the presence of an entropy floor, gener- 
ated by the UV heating, reduces gas densities, and makes 
the H2 molecules in collapsing halos more vulnerable to a 
LW background. Motivated by the above, in this section, 
we study the impact of a LW background on our results. 
We start by a brief discussion of our results without a LW 
background (S I4.1J I. use these results to build up some ex- 
pectations for the impact of the LW background (S I4.2JI . 
and then present the results of simulation runs with LW 
backgrounds (§ 14.31) . 

4.1. Discussion of Results without the LW Background 



13 



As already discussed above, the UV heating produces 
two prominent effects: it boosts the H2 fraction and it 
decreases the gas density. In the case of the individual 
halo studied in Figure [3 and described in § 13.11 and also 
in the case of the analogous halo in the early heating run, 
described in § 13.51 above, we have seen that the overall 
impact of the UV heating is a delay in the development 
of cold dense gas; this delay can be understood by the 
increase in the cooling time, given by the product of the 
two effects above. 

In order to understand the net effect of the UV heat- 
ing on the overall halo population, in Figure^l we show 
the ratio of the H2 fractions (red filled squares), of the 
temperatures (blue empty triangle) and of the mean gas 
density within the central 15pc (red filled circles) in all 
four of our runs with UV heating. Each quantity is com- 
puted for every halo present at z = 23 (or at z = 30 
in the early heating run), in the runs with UV heating, 
and the ratio refers to this value, divided by the same 
quantity in the run without a UV background. 

The figure clearly shows that the H2 fractions are en- 
hanced in a similar fashion in all of the runs (by factors 
ranging from several hundred at low mass, to a few at 
high mass). This is indeed expected: while the H2 abun- 
dance is nearly independent of halo mass in runs with a 
UVB, it is a strongly increasing function of mass in the 
NoUVB case (see Figure 3). Because the "freeze-out" 
value of the H 2 abundance, /h 2 ~2x 1CU 3 is insensitive 
to the background flu x or duration, or to the gas density 
IjOh fc Haimanll2002(K the enhancement factor over the 
NoUVB run is always the same. 

Contrary to the "universal" effect on the H2 fraction, 
the impact of the UV background on the gas density 
depends strongly on the nature of the heating. Not sur- 
prisingly, the flash heating case shows the weakest gas 
dilution; heating the gas for an extended period, at in- 
creasing flux levels, causes larger dilutions. Note that 
the impact on the density tends to diminish for more 
massive halos. This is partly because a fixed amount 
of heating/energy input corresponds to a smaller frac- 
tion of the halo's total binding energy. In addition, the 
H2-cooling time is shorter than 10 7 years in halos with 
M > 10 5 - 8 M© and the UV-heated gas is able to cool 
prior to z = 23. This latter effect is also directly evident 
in the gas temperature ratios, which decrease towards 
larger halos (and decrease below unity). 

The above trends account for the basic results shown 
in Figure El Note that this figure shows only those ha- 
los that develop cold dense gas; i.e. those with M > 
10 5 - 5 M . In the Flash-ionization case, the H 2 fraction 
enhancement in these halos dominates, and results in a 
positive overall feedback. In the Heat0.08 case, the ef- 
fects on the H2 fraction and on the gas density nearly 
cancel each other and the net result is that the UV heat- 
ing has almost no impact. In the Heat0.8 case, the dilu- 
tion of the gas density dominates, and results in a delay 
in the cooling time, and in the development of the cold 
dense gas, by a factor of 1 — 10. 

4.2. The Impact of a Lyman Werner Background - 
Expectations 

The above trends suggest that the UV heating can ren- 
der the halos more susceptible to the negative effect of 
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Fig. 13. — Upper panel: The ratio of the H2 fractions, temper- 
atures, and average gas density between the He at0. 8 and NoUVB 
runs, at 2 = 24 (following the notation in Figure [T2*l , Lower panel: 
The critical value of the background LW flux, Jlw (i n units of 
10 — 21 ergs s — 1 cm -2 Hz -1 sr — 1 ) such that the gas temperature 
cools to 300K by redshift 2 = 18. Only those halos whose gas 
manages to cool by 2 = 18 are shown. Values for the NoUVB run 
are represented by filled black triangles; values for the HeatO.8 run 
are represented by filled red circles. Empty symbols denote the 
critical value of Jlw; such that the gas temperature decreases by 
half between redshifts 2 = 24 and 2 = 18 (see discussion in text). 



a LW background. As argued in lOh k. Haimanl (J2003), 
the H2 photodissociation rate depends linearly on the gas 
density, while the rate of H 2 -forming two-body collisions 
scales with the square of the density; hence density di- 
lution makes H2 photodissociation comparatively more 
important. 

In order to investigate the impact of a LW background 
on the amount of cold dense gas, we performed a set of 
six additional simulation runs. Before describing these 
runs, however, we use the no-LW runs with UV heating 
(Heat0.8) and without heating to develop some expecta- 
tions. These are shown in Figures IT3TI15I as follows. 

In Figure 1131 in the upper panel, we show the ratio 
of the H2 fractions (red filled squares), of the tempera- 
tures (blue empty triangles) and of the mean gas density 
within the central 15pc (red filled circles). The ratios are 
computed in the Heat0.8 and NoUVB runs, as in Fig- 
urc^l but we here use z — 24, rather than z = 23. This 
choice is made to allow for some Compton cooling, but 
to minimize the H2-cooling that occurs after the heating 
is turned off (the latter may not occur if a LW back- 
ground is always on) . In Figure 1141 we explicitly show 
the H2-cooling and Compton-cooling times for each halo 
in the Heat0.8 and NoUVB runs, at z = 24. Note that 
the H2-cooling time is shorter than the Hubble time in 
halos with M > 10 58 M . 

In the bottom panel of Figure 1131 we compute 
the coupled chemical and thermal evolution at fixed 
density, and compute, for each halo, the critical 
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Fig. 14. — Upper panel: H2— cooling time for each halo in the 
Heat0.8 (red, filled circles) and NoUVB (black, empty triangles) 
runs, at redshift z = 24. Lower panel: Compton-cooling timescale 
for the same halos. Note that the H2-cooling time is shorter than 
the Hubble time in halos with M > 10 5 ' 5 M . 



value of the background LW flux, Jlw (in units of 
10~ 21 ergs s^ 1 cm~ 2 Hz -1 sr _1 ) such that the gas tem- 
perature cools to 300K by redshift z = 18. This will rep- 
resent a proxy for a critical value for the LW background, 
above which the halo is prevented from developing cold 
dense gas in the simulation prior to z — 18. The choice 
of the temperature, 300 K, matters relatively little for 
the low-mass halos. On the other hand, we find that the 
critical Jlw we derive for the larger (M > 10 6 M ) halos 
is more sensitive to this choice; in particular, these halos 
have high (> 1000K) virial temperatures, and typically 
never cool down to 300K, even in the absence of a LW 
background. Hence, for these halos, we show (in empty 
symbols) , the critical value of the background LW flux, 
Jlw, such that the gas temperature decreases by half be- 
tween redshifts z — 24 and z = 18. The bottom panel 
in Figure EH shows that the critical Jlw is between 10~ 4 
and 10°, with a relatively large scatter at fixed halo mass. 
There is, nevertheless, a clear impact of the heating by 
the UV background, which reduces the critical LW flux 
by about an order of magnitude (shown as a vertical 
offset between circles and triangles). Note that the ma- 
jority of halos (especially at low masses) never develop 
cold dense gas, and are not shown in the bottom panel. 
In Figure^! we show the ratio of the critical LW back- 
ground as a function of the ratio of the Ff2-cooling time 
(which scales approximately as t coo \ oc T/[n g fn2\) at 
z — 24. Note that there were only 12 halos for which the 
critical LW background was finite in both the Heat0.8 
and NoUVB runs (this excludes the majority of halos, 
which do not form cold dense gas even if Jlw = 0, and 
also those handful of halos that have already formed cold 
dense gas prior to z — 24, in either run). As a result, 
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Fig. 15. — The figure shows the ratio of the amplitude of critical 
LW background (defined to prevent gas cooling between z = 24 
and 2 = 18) between the Heat0.8 and NoUVB runs, as a func- 
tion of the ratio of the H2— cooling times. The scaling is close to 
the Jlw °c t~ ool expected from equating the H2— cooling and H2- 
photodissociation timescales. 



the range shown by this plot is not necessarily represen- 
tative. Nevertheless, the figure shows a clear trend: the 
critical LW background scales nearly as the inverse of the 
cooling time. This can be understood easily: in order to 
prevent the gas from cooling, the H 2 -dissociation time, 
^dissoc ~ 2 x 10 7 (Jlw/10 -3 ) yr, must be comparable or 
shorter than the H2-cooling time. The critical LW flux 
falls somewhat below the value predicted by this scaling, 
because at higher LW fluxes, the H2 abundance starts 
saturating as it approaches its equilibrium value (rather 
than decreasing linearly with time under the influence of 
the background). 

4.3. The Impact of a Lyman Werner Background - 
Simulation Results 

To better quantify the feedback effects of UV heat- 
ing combined with a persistent LW background, we per- 
formed six additional simulations, in which the LW back- 
ground was left on after the UVB was turned off (at 
z — 24.62). Our background flux is constant throughout 
the narrow LW frequency band (11. 18-13. 6eV). We nor- 
malize the specific intensity at the mean photon energy of 
12.87 eV, in units of 10 _21 ergs s _1 cm~ 2 Hz -1 sr _1 . We 
include three values of the LW background: Jlw = 0.001, 
0.01, and 0.1. Each of these three LW backgrounds is ap- 
plied to both our NoUVB and Heat0.8 runs at z = 24.62, 
and is subsequently left on. 

In Figure 1161 we show the resulting CD gas suppres- 
sion, as defined in eq. © and Q). Empty symbols 
refer to runs with no heating (NoUVB), while filled 
symbols refer to runs with heating (Heat0.8). In both 
cases, squares, circles, and triangles denote simulation 
runs with increasing LW backgrounds (Jlw — 0.001, 
0.01, 0.1, respectively). Note that we obtain values of 
o~N,cd{z) < —1 in Figure El this is due to the fact that 
the CD gas disappears from some of the low mass ha- 
los in the presence of strong LW fluxes (c.f the normal- 
ization of eq. @, which tracks relative changes since 

Z = ZUVB,on = 25). 

The results of the simulation runs in Figure 1161 agree 
fairly well with the semi-analytical arguments in § 14.21 
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Fig. 16. — This figure shows the suppression of cold dense 
gas in halos in simulations runs that include a persistent LW 
background. The LW background had specific intensities of 
J LW = 0.001,0.01, or 0.1 (normalized at 12.87 eV, in units of 
10 _21 ergs s — 1 cm -2 Hz -1 sr — 1 ). Each of these three LW back- 
grounds is applied to both our NoUVB and Heat0.8 runs at 
z = 24.62, and is subsequently left on. Values of &m d( z ) (top 
panel) and <5jv c d( 2 ) (bottom panel) are show, as defined in equa- 
tions J7J and I5T 



above. Namely, the value of the LW background at which 
significant suppression occurs by z = 18 in the NoUVB 
runs is found to be approximately Jlw ~ 0.01. In the 
Heat0.8 run, there is significant suppression at z = 18 
already for Jlw ~ 0.001. While the bulk of this suppres- 
sion is due to the UV heating alone (and not the LW 
background; cf. Figure EJ), the LW background does pre- 
vent three additional halos from cooling their gas prior 
to z = 18. This is consistent with the expectation that 
the UV heating lowers the value of Jlw required for ap- 
preciable negative feedback, by a factor of ~ 10. 

More generally, our results reveal that for Jlw < 0.01, 
negative feedback is dominated by UV heating, while for 
Jlw > 0.01, negative feedback is dominated by the LW 
background. Near the threshold value of Jlw ~ 0.01, 
negative feedback transitions from being UV heating 
dominated (< lOOMyr after ZuvB.off) to being LW back- 
ground dominated (> lOOMyr after ZuvB.off)- This 
"transition" behavior can be understood as a combined 
result of two effects: the UV heating is turned off, and its 
impact is transient, as discussed above, while the critical 
LW bac kground scales roughly with the inverse o f the 
density l|Haiman et all 120001 lOh fc Haimanl l200.1l) and 
hence a fixed LW background will have a larger impact 
at lower densities or decreasing redshifts. 

5. CONCLUSIONS 

We used three-dimensional hydrodynamic simulations 
to investigate the effects of a transient ultraviolet (UV) 
flux on the collapse and cooling of pregalactic clouds, 



with masses in the range 10 5 - 10 7 M©, at high red- 
shifts (z > 18). Although in the scenario we envision, 
the radiation is due to nearby PopIII star formation, in 
order to study its effect in a statistical way, we adopted 
a spatially constant but short-lived photo-ionizing back- 
ground throughout the simulation box. This was done 
to mimic the effect of a ~ 100 solar mass star forming at 
z = 25 and shining for 3 Myr. Of course, in reality, the 
closest star can be located at a range of distances and so 
wc effectively covered this range by varying the strength 
of the background. The effect of the ionizing background 
will be strongest on relatively low density gas which is 
in the process of assembling to form halos at later times. 
The sign of the effect has been uncertain with sugges- 
tions of positive feedback due to enhanced H2 formation, 
and negative feedback due to the increased entropy of 
gas in the relic HII region. In addition, we studied the 
combined effects of this transient UV flux and a persis- 
tent Lyman- Werner (LW) background (at photon ener- 
gies below 13.6eV) from distant sources. 

In the absence of a LW background, we find 
that a critical specific intensity of Juv ~ 0.1 x 
10~ 21 crgs s _1 cm~ 2 Hz -1 sr _1 demarcates the transi- 
tion from net negative to positive feedback for the halo 
population. A weaker UV flux stimulates subsequent star 
formation inside the fossil HII regions, by enhancing the 
H2 molecule abundance. A stronger UV flux significantly 
delays star-formation by reducing the gas density, and 
increasing the cooling time at the centers of collapsing 
halos. At a fixed Jtjv, the sign of the feedback also de- 
pends strongly on the density of the gas at the time of 
UV illumination. In either case, we find that once the 
UV flux is turned off, its impact starts to diminish after 
~ 30% of the Hubble time. 

In the more realistic case when a LW background 
is present (in addition to the ionizing source), with 
^lw > 0.01xl0 _21 ergs s^ 1 cm" 2 Hz -1 sr _1 , strong sup- 
pression persists down to the lowest redshift (z = 18) in 
our simulations. Finally, we find evidence that heating 
and photoevaporation by the transient UV flux renders 
the ~ 10 6 M© halos inside fossil HII regions more vul- 
nerable to subsequent H2 photo-dissociation by a LW 
background. 

The results of this study show that the combined neg- 
ative feedback of a transient UV and a persistent LW 
background is effective at high redshift in suppressing 
star-formation in the low-mass halos; this suppression 
will help in delaying the reionization epoch to z = 6 — 10 
as inferred from SDSS quasar spectra and from CMB po- 
larization anisotropy measurements in the 3-yr WMAP 
data. 
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APPENDIX 

COMPARING SEMI-ANALYTIC AND SIMULATION MASS FUNCTIONS 

Although it is not directly relevant to the radiative feedback processes analyzed in the body of this paper, it is 
interesting to examine the mass function of dark matter halos we find in our simulations, and compare these to semi- 
analytic models. Such a comparison is especially interesting, since our simulation corresponds to a biased region that 
is overdense on the scale of the simulation box (rather than fixed to the cosmic mean density) . We are not aware of a 
previous study of the halo mass function derived in such a biased region. Here we present a preliminary comparison, 
and postpone more detailed work for a future paper. 

According to extended Press-Schechter formalism (EPS), the contribution of halos with masses greater than M m i n 
to the mass fracti on inside regions of mass scale M=( 4/3)7ri? 3 (/0n) and extrapolated 3 mean overdensity Sq, can be 
expressed as (e.g.. lBond et al.lll99lULacev k. Coleill993j) : 



^col(> M„ 



= erfc 



S c (z) - S 



^/2[S(M min ) - S(M )] 



(Al) 



where S c (z) is the critical linear overdensity at halo virialization, and S(M) is the variance of the present-day linear 
overdensity field filtered on scale M. For comparison purp oses below, we chose M m ; n = 10 5 M Q . 

To compare our numerical results with equation I)A1|I . we chose our biasing scale, Mq, such that the enclosed 
spherical volume is equal to the volume of the refined region in the simulation, (4/3)-7ri?g = Z 3 of j ncd = 1/64 (h^ 1 Mpc) 3 , 
or equivalently, a mass scale of Mq = 6.4 x 10 8 Mq. Keeping in the spirit of the linear, Lagrangian nature of EPS, we 
obtain our overdensity bias, <5o> from the initial redshift of our simulation, z- ln i t — 99, before the refined region in our 
Eulerian code has a chance to become too contaminated by inflow from outside particles. We obtain S(zi n it) — 0.1637 
in our refined region, making our linear overdensity bias Sq = S(z- m i t )/ D(z[ n it) — 12.8. 

We sho w th e comparison of this biased mass fraction with our simulation in Figure lATl The solid line is determined by 
equation IjAljI . using the values quoted above. The points denote the values obtained from the cosmological simulation 



3 We adopt the standard convention to work with the density field linearly extrapolated to z - 
is the linear growth factor normalized so that D(z = 0) = 1 (e.g., Liddle et al. 1996). 



: 0, i.e. S(z = 0) = 8(z)/D(z), where D(z) 
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Fig. Al. — Fraction of mass of the biased refined volume con tain ed in halos with masses greater than 10 5 Mq. The solid line is determined 
from the extended Press— Schechter formalism, using equation 1AH , The points denote the values obtained from the cosmological simulation 
by su mming over halos found with the HOP algorithm, in the case of no background radiation field. The dashed line is obtained by equation 
IA1I . but for a mean sample of the universe (i.e. <5q = 0, S(Mq) = 0). 



by summing over halos found with the HOP algorithm, in the case of no background radiation field. The dashed line 
is obtained by equation IjAlll . but for a mean sample of the universe (i.e. <5 = 0, S(M n ) = 0). 

We note that mass fractions obtained f rom the simulation are a factor of ~ 2 lower than those obtained by EPS. To 
further probe this discrepancy, in Figure lA"2l we plot the analogous mass functions at z = 25 (top panel) and z = 20 
(bottom panel). It is evident from the figure that abundance of low-mass halos is over-predicted by EPS with respect 
to the simulation, while the abundance of the largest halos in the simulation fits the EPS prediction fairly well. 

This discrepancy might be due to several reasons. Firstly, it is already known that EPS mass-functions in the 
low-redshift regime suffer from similar over-predictions of low-mass halos, as well as an under-prediction of high- 
mass halos, where "low" and "high" are defined with respect to the characteristic collapse scale at z (Jenkins ct al. 
2001: ISheth fc Tormenl 119991) . with the mass functions differing by up to a factor of ~ 1.6. Another contributing 
factor to the discrepancy could be the fact that the refined region of our simulation is a cubical perturbation, while 
parameters in the standard EPS are derived assuming spherical perturbations of the real-space density field. Finally, 
it has been shown that high-resolution cosmological simulations are too small to provide accurate mass fu nctions at 
high redsh ifts, since they ar tificially cut-off the density modes larger than their box sizes ( Barkana & Loeb 2004) . In 
particular. iBarkana fc Loea IJ2004J) show that at z = 20, the true cosmic mean mass function can be a factor of several 
higher than would be derived from a (1 Mpc) 3 simulation box, with periodic boundary conditions normalized to the 
mean density. Assuming that our biased simulation box suffers from a similar un derestimate, the mass fu nction could 
b e consistent wi t h the EPS prediction (after applying the correction proposed bv ISheth fc Tormenll 1999ft . 

Yoshida et alJ l)2003f) ha ve obtained a good fit a t high redshifts (z ~ 20) between EPS mass functions and those 
obtained from simulations. IBarkana fc Loebl (2004) show that after correcting for the above missing large-scale power, 
that their results are consistent with the EPS mass function with the Sheth-Tormen correction. However, these results 
describe an unbiased simulation at mean density. We present the first comparison of the biased EPS and numerical 
mass functions at high redshifts. A detailed study on such a comparison is beyond the scope of this paper, nor does 
it have an impact on our main results. 
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Fig. A2. — Mass functions for the models analogous to Figure lAll at z = 25 (top panel), and z = 20 (bottom panel). 



